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Abstract: In Q it was shown how non-perturbative corrections to gravity can resolve the 
big bang singularity, leading to a bouncing universe. Depending on the scale of the non- 
perturbative corrections, the temperature at the bounce may be close to or higher than the 
Hagedorn temperature. If matter is made up of strings, then massive string states will be 
excited near the bounce, and the bounce will occur inside (or at the onset of) the Hagedorn 
phase for string matter. As we discuss in this paper, in this case cosmological fluctuations 
can be generated via the string gas mechanism recently proposed in 0. In fact, the model 
discussed here demonstrates explicitly that it is possible to realize the assumptions made in 
Q in the context of a concrete set of dynamical background equations. We also calculate 
the spectral tilt of thermodynamic stringy fluctuations generated in the Hagedorn regime 
in this bouncing universe scenario. Generally we find a scale-invariant spectrum with a red 
tilt which is very small but does not vanish. 
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1. Introduction 

Recently, a new structure formation scenario has been put forward in |^, [| (see also || 
for more in depth discussions). In these works it was shown that string thermodynamic 
fluctuations in a quasi-static primordial Hagedorn phase ^] ( during which the temperature 
hovers near its limiting value in string theory, namely the Hagedorn temperature, Th, fflj, 
the metric is almost static, and hence the Hubble radius is almost infinite) can lead to 
a scale-invariant spectrum of metric fluctuations. To obtain this result, several criteria 
for the background cosmology need to be satisfied (see ||]). First of all, the background 
equations must indeed admit a quasi-static (loitering) solution. Next, our three large 
spatial dimensions must be compact. It is under this condition that |ll[] the heat capacity 
Cy as a function of radius r scales as r 2 . Thirdly, thermal equilibrium must be present over 
a scale larger than 1mm 1 during the stage of the early universe when the fluctuations are 

1 This scale is obtained by evolving our present Hubble radius back to temperatures of the order of the 
scale of Grand Unification using the equations of Standard Cosmology. 
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generated. Since the scale of thermal equilibrium is bounded from above by the Hubble 
radius, it follows that in order to have thermal equilibrium on the required scale, the 
background cosmology should have a quasi-static phase. Finally, the dilaton velocity needs 
to be negligible during the time interval when fluctuations are generated. 

It is not easy to satisfy all of the conditions required for the mechanism proposed in 
to work. In the context of a dilaton gravity background, the dynamics of the dilaton is 
important. If the dilaton has not obtained a large mass and a fixed vacuum expectation 
value (VEV) at a high scale, then it will be rolling towards weak coupling at early times. 
This will lead to 0, a phase in which the string frame metric is static, and thus the 
string frame Hubble radius will tend to infinity, i.e. H ~ 0. However, the large dilaton 
velocity during this phase changes the spectrum of fluctuations coming from the stringy 
Hagedorn phase || g|. In addition, the duration of the Hagedorn phase is too small for the 
establishment of thermal equilibrium on sufficiently large scales || (also see [10 j). These 
problems can be overcome in a model which has fixed dilaton and admits a bouncing solution 
2 . Such a bouncing solution for the Einstein frame metric turns out to be a key prediction 
of a higher-derivative theory of gravity recently proposed by a subset of the current authors 
IQj. Here, we will show that in the model proposed in jlj], the conditions required for the 
structure formation mechanism of [^ to work are naturally satisfied. 

Obtaining a bouncing cosmology without introducing any pathologies is a challenging 
task. It is believed that superstring theory will lead to a resolution of this problem by 
modifying the ultraviolet behavior of the Einstein-Hilbert action. Any generally-covariant 
modification of gravity involves higher derivative extensions of the Einstein-Hilbert action. 
Thus it is no surprise that string theory indicates such corrections [14, 15]. However, higher 
derivative theories usually contain ghosts, thereby making them inconsistent. There are 
two different mechanisms known which can avoid this problem. The first method exploits 
the fact that gravity is a gauge theory and therefore, even if there are ghosts, they may be 
benign [pi]]. It then becomes possible to avoid the "dangerous ghosts" by choosing special 
combinations of the higher derivative corrections. Gauss-Bonnet gravity |15| is a familiar 
example of this class (but see also [[17]] for a more critical analysis) . A second method relies 
on non-perturbative physics and is equally applicable to gauge theories (like gravity) or 
non-gauge theories, such as scalar field theories. The way non-perturbative modifications 
avoid the problem of ghosts is by ensuring that the corrections do not introduce any new 
states, ghosts or otherwise. (Perturbative corrections, i.e. up to some finite order in higher 
derivatives, inevitably introduce new states.) A well known example of such a theory is 
the p-adic string theory where the kinetic part of the action contains an infinite series of 
higher derivative terms |l8| (see also | IS , [2(| ) . 

A non-perturbative ghost-free model of gravity containing "string-inspired" higher 
derivative corrections was recently proposed in Q, it was also shown that such theories 
offer a possible resolution to the big bang singularity as they admit bouncing solutions. 
The discussion of (lj was in the context of radiative matter. However, typically the energy 



2 A bouncing cosmology automatically implies a phase of acceleration. However, in general this phase 
will be short (of the order or smaller than one Hubble expansion time), and will thus not lead to a long 
period of cosmological inflation. 
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density involved during the bounce is of the order of the cut-off scale which could be of the 
order of the string scale. It is then natural to imagine that for a string-scale energy density 
one would excite at least the fundamental strings. If we describe the bounce in terms of 
three spatial dimensions, assuming that the other spatial dimensions are stabilized, then 
one could naively expect a gas of closed strings in a Hagedorn phase during the bouncing 
phase. Here, we first generalize the analysis of [|IJ and show that non-perturbative gravity 
actions admit bouncing solutions and therefore can resolve the big bang singularity for 
general forms of matter (not just radiation), a stringy Hagedorn phase for instance. Then, 
we show how the structure formation scenario of ||] can be implemented. 

As we shall argue, to generate the appropriate metric fluctuations, we would require 
the bounce to occur just after the onset of the Hagedorn phase, when although the fluctu- 
ations in energy is dominated by closed string excitations, the energy density itself is still 
dominated by the massless string modes (radiation) 3 . There are three reasons. First, in 
order to match the current observations, the amplitude of the metric fluctuations, which 
will be determined by the ambient temperature, must be small. It turns out that this is 
easier to achieve if the bounce occurs close to the onset of the Hagedorn phase. Secondly, 
if the energy density is dominated by the Hagedorn gas of closed strings, then it is hard 
to keep the dilaton fixed, because in the Einstein frame strings couple to the dilaton and 
will therefore source it. Finally, in order to apply local physics we need to define a local 
density which can only be done if energy goes as volume, as it does for radiative matter. 
The crucial point to note is that all the conditions required for the success of the model 
are met near the onset of the Hagedorn phase. 

The paper is organized as follows: In section § we first review non-perturbative gravity 
models and then discuss the background bouncing geometries that emerges in these mod- 
els. In section || we discuss the spectrum of stringy perturbations that can be generated 
in the bouncing phase if the bounce occurs within the Hagedorn phase. In particular we 
estimate the spectral tilt and amplitude for some specific cases. We conclude in section [| 
by summarizing the new structure formation mechanism based on singularity free bounc- 
ing universe scenario, and also point out possible caveats in the model. There are three 
appendices devoted to discussing, the technical details of obtaining the bouncing solutions 
(appendix |A|), the thermodynamics near the onset of the Hagedorn phase (appendix [B]), 
and the dynamics of possible light moduli fields during the course of the bounce (appendix 

0). 

2. Non-Perturbative Gravity and Bouncing Universes 
2.1 Review 

In string theory, higher-derivative corrections to the Einstein-Hilbert action appear already 

3 Briefly (see appendix |b|) , the reason this happens is because the energy fluctuations are proportional 
to the heat capacity which in the Hagedorn phase goes as Th/(Th — T), while the energy only goes as a 
logarithm of this ratio. Therefore, as the temperature reaches close to the limiting Hagedorn temperature 
the heat capacity corresponding to the excitations of the massive modes catches up with those of the 
massless modes much faster than the energy itself. 
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classically (i.e. , at the tree level), but we do not preclude theories where such corrections 
(or strings themselves) appear at the loop level or even non perturbatively. From string 
field theory (either light-cone or covariant) the form of the higher-derivative modi- 
fications can be seen to be Gaussian, i.e. there are e a factors appearing in all vertices 
(e.g., (e'-'c/)) 3 ). These modifications can be moved to kinetic terms by field redefinitions 
(4> —* e~ (f>). The nonperturbative gravity actions that we consider here will be inspired 
by such stringy kinetic terms In this context we note that similar non-local modifica- 
tions of gravity has been considered previously in the literature |2^| , while only recently 
possible modifications to the evolution during the Hagedorn phase due to (perturbative) 



higher derivative corrections was considered in [23]. 

It was also noticed in Ref Q that if we wish to have both a ghost free and an asymptot- 
ically free theory of gravity , one has little choice but to look into gravity actions that are 
non-polynomial in derivatives, such as the ones suggested by string theory. Thus we start 
by considering the simplest non-polynomial generalization of the Einstein-Hilbert action 

S = Jd 4 x ^F(R) , (2.1) 

where 

oo 

F(R)=R+j2ijk RnnR > ( 2 - 2 ) 

n=0 * 

and M* is the scale at which non-perturbative physics becomes important. c n 's are typically 
assumed to be ~ 0(1) coefficients. It is convenient to define a function, 

/ 00 / \ \2(i+l)\ 
!■(*)= (l-.J>(£) ). (2.3) 

One can roughly think of p 2 T(— p 2 ) as the modified inverse propagator for gravity (see |i 
for details). In Q it was shown that if T(A 2 ) does not have any zeroes, then the action 
is ghost-free 5 , and in fact one introduces no other states other than the usual massless 
graviton (there are no extra poles in the propagator). One can also check that if T goes 



4 While perturbative unitarity requires the theory to be ghost free, in order to be able to address the 
singularity problem in General Relativity, it may be desirable to make gravity weak at short distances, 
perhaps even asymptotically free p3{. 

5 One may also worry about classical instabilities which usually plague higher derivative theories, gener- 
ically they go by the name of Ostrogradski instabilities (see |2f| for a review). However, these are basically 
the classical manifestations of having ghosts in your theory (non-unitarity or having negative norm states 
can be exchanged for having arbitrarily large negative energy states, which classically manifests EIS El CcltcLS - 
trophic instability). Since the non-local theories under consideration do not contain any ghosts we also do 
not expect to find such instabilities. One way to see how these theories may avoid the Ostrogradski insta- 
bility argument, which is valid for finite higher derivative theories, is by noting that one cannot construct 
the usual Ostrogradski Hamiltonian because one cannot identify a "highest derivative" in such non-local 
actions. Also, in arriving at the Ostrogradski Hamiltonian, one assumes that all the derivatives of the 
field (except the maximal one) are independent canonical coordinates. This is no longer true for theories 
with infinite derivatives. For instance, one cannot independently choose an initial condition where one can 
specify arbitrarily the values of all the derivatives of the field. See ]l8[ for a nice discussion regarding this. 
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like an exponential, as expected from stringy arguements, then not only is it ghost-free, 
but it also describes and asymptotically free theory. 

Since we are interested in cosmological solutions, and in particular homogeneous and 
isotropic cosmologies with a metric (for simplicity we consider the case of a spatially flat 
universe) given by 

ds 2 = -dt 2 + a{t) 2 (dx 2 + dy 2 + dz 2 ) , (2.4) 

where a(t) is the scale factor, it is sufficient to look at the analogue of the Hubble equation 
for the modified action (2.1,2.2) 6 . This equation is |]26| , 



~ F 1 °° 3 °° 

Goo = F R 00 + --Fo- 00 -aF --J2 ^™ - - J] F n (a n ^R) = T 00 , (2.5) 

n=l n=l 

where we have defined 

F m = V u n ~ m ^- . (2.6) 

Q\Jn R V J 

n=m 

It was shown in that ( |2.5[ ) admits exact bouncing solutions of the form 

a(t) = cosh . (2.7) 

in the presence of radiative matter sources and a non-zero cosmological constant. Note 
that this is an exact result, but other approximate bouncing solutions can now also be 
constructed in the presence of arbitrary matter sources, which can match the gradual 
contraction before and expansion after the bounce. This is important for several reasons. 
Firstly, it shows that the cosmological constant is not necessary for the bounce as was 
already argued in [|J . Secondly, it shows that one can have bounce not only in the presence 
of radiation but also other sources of matter, such as a stringy hagedorn phase which is 
expected in the early universe. Finally, we can obtain bouncing solutions for a much larger 
class of non-perturbative actions thereby exhibiting the intrinsic quality of these corrections 
to resolve singularities. We now discuss these approximate solutions. 

2.2 Approximate Bouncing Solutions for Generic Matter Sources:Energy and 
Distance Scales 

For the purpose of this paper, since we are interested mainly in the bouncing phase, we will 
assume the bounce to be of the form Q2.7D without loss of any generality. The crucial point 
is to note that the entire dynamics can be analyzed by separating the time evolution into 
three different regimes. Very close to the bounce all the non-perturbative terms become 
important and one can approximately compute both the left and the right hand side of 
the modified Hubble equation (|2.5| ) to obtain solutions which are valid close to the bounce, 
i.e. for A|t| < 1. However, it is easy to check (see Appendix A for details) that the non- 
perturbative corrections fall faster than the Einstein tensor as one goes away from the 



6 Just as in ordinary Einstein gravity, here also the Bianchi identities (conservation equation) ensure that 
for FRW metrics the Ga equation is automatically satisfied when the Goo equation is. 
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bounce, and that around A|i| ~ 1, Einstein gravity takes over. Thus for X\t\ > 1, both in 
the contracting and expanding phase, we have the usual evolution governed by ordinary 
gravity coupled to matter. The main task therefore is to understand the evolution near 
the bounce (if there exists one). 

For this purpose, we will not assume any specific form for the energy density, namely 
the source term, see the right hand side of (2.5), but just that it is a function of the scale 



factor: 

Too = p(a) = p(a(t)) . (2.8) 

What we now can do is to expand both the left and right hand side of ( |2.5[ ), i.e. , Goo 
and p(a(t)) as a power series in t 2 for the bouncing ansatz (2.7). Then by matching the 



coefficients of the first few powers in t 2 , we can determine the dynamical parameters (see 
appendix [A| for details). In principle this calculation can be done for any matter density, 
as long as we know p(a), but for the cases of practical interest p{a) behaves as an ideal gas 
with a specific equation of state parameter u> 

P = uj P , (2.9) 

where p is the pressure of the fluid. This implies 

p = Ph a~ Z{1+w) (2.10) 

where pb is the maximal density attained at the bounce point. Assuming such an ideal 
gas approximation, after some algebraic manipulations one finds that the time scale of the 
bounce, A -1 , is determined by an equation of the form 

K (l)^° (2 ' u) 

where K, is a dimensionless function determined entirely in terms of T (see appendix |A|) i.e. 
, the non-perturbative corrections. This is not surprising as it is these corrections which 
are responsible for causing the bounce in the first place. Also, one straight away sees that 
since the coefficients in the function K, are typically expected to be ~ 0(1), we should have 

A ~ M* . (2.12) 

unless there is some fine-tuning involved. Next, we find that the maximal energy density, 
Pb , is given by an equation of the form 

3M 2 \ 2 G(^) 

Pb= g ' >0- (2.13) 

where Q is again a dimensionless function which can be expressed in terms of T, see appendix 
|A| for details. It's specific expression is complicated and not very illuminating, but again 
one does not need to know the details to understand the basic physics. For instance what 
is important is to be able to determine the scale of energy density at the bounce point. 
Again, if no fine-tuning is involved, then Q ~ 0(1) and we have 

p ~ M 2 X 2 ~ M 2 M 2 . (2.14) 
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It is now clear when we can start to see some effects of a stringy Hagedorn phase: 
During the phase of radiation domination, the density is given by p ~ T 4 . Therefore we 
would expect stringy effects to show up if pb > T]j ~ Mf, or in other words as long as 

M 2 

M, > j±. (2.15) 

If the higher derivative corrections are genuinely stringy in origin, then M* ~ M s , and one 
sees that ( |2.15| ) is easily satisfied as the string scale is expected to be less than the Planck 
scale. Thus, in general, one would indeed expect our bouncing universe to probe the very 
stringy Hagedorn phase of matter. 

However, as we will see later, for phenomenological reasons we will be interested in the 
regime where the bounce probes just the onset of the Hagedorn phase. This is possible if 
one allows for some fine-tuning in either ( |2.11| ) or ( p,13|) . For instance if A <§C M* ~ M s , 
then one can easily have pb ~ Tfj ~ Mf as the last "~" in ( p. 14 ) is no longer true. 



Finally, we find that in order to have a consistent bouncing solution we also have a 
constraint equation 7 which looks like 

1 + u) = > (2.16) 

(J- is again a dimensionless function expressible in terms of T). For a given equation 
of state, together with ( |2.11| ) it over-determine A 2 , thereby restricting the set of non- 
perturbative actions (T's) which can exhibit bouncing solutions of the nature (|2.7f), for 
a given equation of state. However firstly, this shows that non-perturbative actions can 
resolve the singularity not only in the case of radiation (plus a cosmological constant), but 
also for other general matter sources. In particular one can check that for a large number 
of ghost free actions one can find solutions with u = 1/4, that is for pure radiation, while 
we can also find actions admitting solutions deep inside a stringy Hagedorn phase where 
energy remains approximately constant, i.e. ui = 1/3. 

Secondly, since we have the freedom to choose uj, even while restricting ourselves to the 
dominant energy condition (the inequalities in ( |2.16| , |2.13j ) which essentially implies that for 
the A which solves the equations, T, Q have to be positive), we can find a much larger class 
of actions which manifest the bouncing behaviour, at least for some uj. This illustrates the 
intrinsic capability of non-perturbative higher derivative corrections to resolve the big bang 
singularity problem. We note, in passing, that for usual Einstein gravity or R 2 gravity, the 
dominant energy condition cannot be satisfied. Indeed, one does not have a non-singular 
bounce in those theories. 

Before concluding this section, let us briefly discuss about another important physical 
scale, that of the size of the universe, Rb, at the bounce point. In the context of a bouncing 
universe, it is quite natural to have this scale much larger than the string scale. In order 
for the structure formation mechanism of H to work, we require that Rb is larger than 



7 This is not surprising, as such constraints were found when dealing with the exact solutions in pj, and 
in fact coincide for these special cases. 
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1mm. In the context of the string gas cosmology setup proposed in || |2|, this is not 
natural. Obtaining such a large initial universe in the context of string gas cosmology 
leads to a similar entropy problem as is faced in Standard Cosmology (for a possible 
solution of this problem making use of a gas of branes, see [27|). An interesting aspect of 
our bouncing universe proposal is that there is no entropy problem, as long as the universe 
in the contracting phase started out large (which, in the context of a cold initial state, is 
quite natural). 

One more comment: so far we have not discussed the role of the dilaton 0. We have 
tacitly assumed that the dilaton has been fixed, such that the string coupling 9s ~ e <$>l M P is 
constant. In fact, we can allow the dilaton to be dynamical near the bounce, see Appendix 
|C^; typically the dilaton attains a constant value before and after the bounce and jumps by 
a finite amount during the course of the bounce. The amount of jump however depends 
on the initial conditions and can be made to be very small, so that it becomes a good 
approximation to treat the dilaton profile near the bounce to be a constant. 



3. Hagedorn Phase and Structure formation 
3.1 The New Cosmological Scenario 

The bouncing solution obtained in the previous section has interesting consequences for 
cosmology if the energy density at the bounce point is large enough to lead to a Hagedorn 
phase of string matter. 

We begin with a sketch of the space-time diagram obtained in our scenario (Fig. 1). 
The vertical axis in Fig. 1 denotes physical time, with the origin of time chosen such 
that t = is the bounce point. The horizontal axis denotes comoving spatial coordinates. 
The two dotted vertical lines labeled by k\ and hi denote the comoving wavelengths of 
cosmological fluctuations with comoving wavenumber k\ and &2- The solid curve is the 
Hubble radius, defined as "H -00 , where Ti is the Hubble expansion rate in conformal time. 
The key point is that the Hubble radius becomes infinite at the bounce point. In particular, 
fluctuations on all scales of interest in cosmology today are sub-Hubble in a time interval 
about the bounce point. The evolution of the Hubble radius as a function of time for t > 
is analogous as in the string gas cosmology setup discussed in Q , where it was shown that 
string thermodynamic fluctuations during an early quasi-static Hagedorn phase can lead 
to a scale-invariant spectrum of curvature fluctuations. 



Let us recall the origin of the Hagedorn phase 0: In the context of perturbative string 
theory, there is a maximal temperature Th, the Hagedorn temperature J7|, for a string 
gas in thermal equilibrium. As the universe contracts in our scenario, the temperature 
initially increases as in the usual radiation phase of standard cosmology. The energy is in 
the radiative degrees of freedom, the momentum modes of the strings. As the temperature 
approaches Th, it becomes possible to excite the string oscillatory and winding mode 
degrees of freedom. In fact, close to the Hagedorn temperature, most of the energy in the 
string gas is contained in winding strings [28, 2£j. 
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Figure 1: Space-time diagram (sketch) showing the evolution of a fixed comoving scale in our 
bouncing cosmology under the assumption that the bounce phase lies entirely within the Hagedorn 
phase. The vertical axis is time, the horizontal axis is comoving distance. The bounce time is 
chosen to be t = 0. The Hagedorn phase corresponds to the time interval — tn < t < tu, the 
bounce phase is the interval — ts < t < tB- Outside of the bounce phase, the Hubble radius evolves 
as in a radiation phase of standard cosmology The black solid curve represents the Hubble radius 
if -1 . The red vertical line labelled by k depicts a fixed comoving scale. This scale is smaller than 
the Hubble radius in the central region of the bounce, exits the Hubble radius at the time t i (k), 
and re-enters at a late time tf(k). 



The crucial point of the mechanism of || is the following: provided that our three 
large spatial dimensions are compact, the heat capacity Cy of a gas of strings in thermal 
equilibrium scales as r 2 with the radius of the box. The heat capacity determines the root 
mean square mass fluctuations via 

((5M) 2 ) = T 2 C V , (3.1) 



which then in turn induce metric perturbations. In longitudinal gauge (see e.g. [30, 31] 



for reviews of the theory of cosmological perturbations), and in the absence of anisotropic 
matter stress at late times, the metric takes the form 

ds 2 = -(1 + 2<S>)dt 2 + a(t) 2 (l - 2$)dx 2 , (3.2) 

where t is physical time, x are the comoving spatial coordinates of the three large spatial 
dimensions, a(t) is the cosmological scale factor and 3>(x, t) represents the fluctuation 
mode. On scales smaller than the Hubble radius, the metric fluctuations are driven by the 
matter fluctuations, on super-Hubble scales it is the metric fluctuations which are dominant 
since the matter oscillations have frozen out (thermal fluctuations cannot keep up with the 
Hubble expansion). Thus, to compute the late-time spectrum of metric fluctuations, the 



philosophy adopted in [2| (see also |32|] for a corresponding treatment of fluctuations in 
inflationary cosmology and |j]J for a discussion on the spectrum of thermal fluctuations 
corresponding to ordinary particles) is, for inhomogeneities on a fixed comoving scale, to 
follow the matter fluctuations until the scale exits the Hubble radius at time U(k), and to 
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use the Einstein constraint equations to determine the metric fluctuations at that time. 
The metric fluctuations will then be conserved on super-Hubble scales. 

As long as the background dilaton velocity is negligible, the Einstein constraint equa- 
tions (specifically, the relativistic Poisson equation) 

V 2 $ = AirGSp, (3.3) 

(where p is the energy density) then imply that the power spectrum of the metric fluctu- 
ation variable is scale-invariant (for details on how to get from the spectrum of mass 
fluctuations to the power spectrum of $ see 0). 



In order for (3.3) to apply, it is important that metric and matter fluctuations are 



related via the Einstein equations, and not e.g. via the dilaton gravity equations [12] 
(see || for a detailed discussion of this point). In our bouncing cosmology, the dilaton 
is taken to be constant throughout the evolution of the universe. Thus, there are no 



dilaton corrections to (|3.3|). We also note that the higher derivative corrections to (3^3) are 
suppressed by (V 2 /M 2 ) ~ (k/M*) 2 (see Q for a quantitative discussion). The comoving 
scales that we are interested in corresponds to the Hubble radius today, and one can easily 
check that this ratio is hugely suppressed (see section [B] for instance) in a non-inflationary 
scenario such as ours. In addition, the fact that there are no additional dynamical degrees 
of freedom compared to those in Einstein's theory ensures that we need not worry about 
extra unstable (ghost like) or spurious excitations of the metric. 

We are thus led to the following structure formation scenario: The universe begins in 
a cold contracting phase. The initial Hubble radius is set by the initial gas density, and it 
is natural to assume that the size of space is at least as large as the initial Hubble radius. 
If we begin the evolution at a temperature lower than the present radiation temperature, 
we are then ensured that the size of space at the bounce is sufficiently large to contain all 
fluctuation modes of current interest in cosmology 8 . Local processes in at the beginning 
or during the course of the contraction phase establish thermal equilibrium over a scale 
which contains all fluctuation modes of interest to us. Matter remains homogeneous and 
isotropic throughout the bounce. In particular, during the Hagedorn phase of the bounce, 
we have a gas of strings. In the same way that a network of strings forming at a phase 



transition in a field theory model with cosmic strings [34] will contain strings of length 
larger than the Hubble radius (in particular, it will contain winding modes), after the 
transition to the Hagedorn phase, our gas of fundamental strings will contain winding 
modes, and it is the presence of these winding modes which are the key ingredient to 
obtain the form of the specific heat we are using. Thus, we are able to use the results of 
string thermodynamics to compute the spectrum of matter perturbations which then seed 
the scalar metric fluctuations. 

3.2 Thermodynamics, a Brief Review 

In order to compute the power spectrum of density fluctuations we have to compute the 
fluctuations in string energy, AE, for these closed string modes inside an arbitrary volume 

8 Note, in particular, that the bounce radius Rt, is not set by the physics which determines the bounce, 
but is in fact much larger than the length scale of the local bounce physics. 
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V. Once we obtain AE(V), there is a well defined prescription on how to go over to 
the power spectrum V{k) as described in the next subsection. In order to obtain AE and 
eventually the power spectrum we have to understand the thermodynamics of the Hagedorn 
phase and here we start with a brief review of the same. 

Specifically we are interested in a situation where we have d = 3 large but compact 
spatial dimensions. The entropy S of a string gas in thermal equilibrium in this background 
was derived in [11] (also see ||, |4|). During the Hagedorn phase, and in the case of a volume 
which is large in string units, the result is 



S = £ + bv 



£ 



2d-l 



(2d - 1) 



exp 



c(£ — kv) 
^273 



where we have introduced the dimensionless variables 



£ = — — and v 



VT, 



H 1 



(3.4) 



(3.5) 



for convenience. V is the volume of the d large dimensions, and b, c, k are 0(1) constants 
(see |2^, H|, |]]). This formula is valid in the regime £ > kv (which is characteristic of the 
Hagedorn phase) and v 3> v 2 ^ (which means that the volume of the universe is much 
larger than the string length). 

Once we have the microcanonical entropy S(E, V) as a function of the energy and the 
volume, it is straight forward to compute all the thermodynamic quantities. In particular 
we find that the inverse temperature (5 is given by 



p dE 



which implies 



0H 



1 + 



c£ 



2d-l 



(2d-l)h; 2 /3 



exp 



c(£ — kv) 
^2/3 



where we have used Using the same approximation we find from ([■ 

specific heat (or rather the heat capacity) is given by 



d£_ 
df 



(2d-l)\v 2 / 3 
cT(T H - T) ' 



(3.6) 

(3.7) 
that the 

(3.8) 



which gives the characteristic scaling Cy ~ r 2 as a function of the radius of the volume 
which, by the discussion in the previous subsection implies a scale-invariant spectrum of 
metric fluctuations. 

If the decrease in the Hubble radius at the end of the Hagedorn phase is not instan- 
taneous, then the time U(k) (and thus also the corresponding temperature T(ti(k)) will 
depend on k. This will lead to a red tilt in the spectrum. In section ( |3.4|) we will see that 



using (3/7) and (3^8) one can calculate the spectral tilt in the bouncing universe scenario 
quite precisely. 
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3.3 Spectrum of Perturbations 

The dimensionless power spectrum of scalar metric fluctuations is denned by 

P*(fc) = k 3 \$(k)\ 2 (3.9) 

where is the k'th Fourier mode of <J?(x), and we are using the convention of defining 
the Fourier transform of a function f(x) including a factor of the square root of the volume 
of space. 



Making use of the relativistic Poisson equation (3.3), (3.9) becomes 



ffe(fc) = 16TT 2 G 2 k- 1 \6p(k)\ 2 . (3.10) 

The Fourier space density perturbation of wavenumber k can be determined from the mean 
square mass fluctuation |<5M| 2 (r) in a sphere of radius r(k) = k^ 1 via 

\5p(k)\ 2 = k 3 \5M\ 2 (r(k)), (3.11) 

and hence ( |3.9| ) becomes 

P*(ife) = 16ir 2 G 2 k 2 (5M) 2 (r(k)) . (3.12) 



Plugging the expression (3.8) for the specific heat capacity of a string gas into (|3.ip and 



then into ( [3.12 ) we finally obtain the following result for the power spectrum ||, [5j]: 



P»(k) = im^c-^T^krik)) 2 ^ , (3.13) 

where for each value of k, the temperature T is to be evaluated at the time that the mode 
k exits the Hubble radius at time ti(k). 

We see immediately that as a consequence of the specific stringy scaling of the heat 
capacity, we obtain a scale-invariant spectrum of fluctuations in the approximation when 
T(ti(k)) is independent of k as it will be if the Hubble radius drops abruptly. If the 
decrease of the scale factor is not instantaneous, then a small red tilt of the spectrum will 
be induced since smaller length scales exit the Hubble radius later, and the temperature 
then is smaller. We will come back to the issue of the magnitude of the spectral tilt shortly. 

The overall amplitude of the spectrum is, as in 0, ||, determined by two factors, firstly 
the factor G 2 T^ which is the fourth power of the ratio of the Planck length and the string 
length, and secondly the final factor in ( 3.1 3| ) which depends on how close the temperature 



at the bounce point is to the Hagedorn temperature. 
3.4 Spectral Tilt: General Formula 

In this subsection we compute the spectral tilt rj s — 1 which is defined by 

P(k) ~ k^- 1 . (3.14) 

To keep things general we do not specify the function T = T(a). The spectral tilt can 
be calculated straightforwardly from ( 3.13| ). The key factor is the factor 1 — T/Th in the 
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denominator of the expression. This is the term that carries the dominant dependence on 
k. The reason is that at all times ti(k), the temperature is very close to the Hagedorn 
temperature so that the only other factor which depends on k, namely the factor T in the 
numerator of ( 3.1 3| ) can be taken to be constant. From ( 3.13j ) 

In 7? = InC -ln(l-T/T H ), (3.15) 

where the constant C stands for the product of all terms in (|3.13| ) which do not depend on 
k. Hence, 

dlnV _ dlnV da 

dk da dk ' 

where the comoving wave- vector k exits the Hubble radius when the scale factor a is given 
by 

k = Ha = a (3.17) 

which implies 

k = a. (3.18) 



Note that ( 3.17 ) determines A; as a function of the scale factor a, and through its dependence 
on time, also as a function of time, k = k{a) = k(a(t)). 
From ( j3~T7D we find 

da l da \ .. f/i,- a t . 

(3.19) 



dk 



( da\ .dk 
\dt) '~dt 



Thus, 



dlnV dlnV a 



dk 



Now, straightforwardly one finds 



dlnV 



da d 
T'T H 



(3.20) 



(3.21) 



da ~ (Th — T)T ' 

where ' denotes the differentiation with respect to the scale factor. Putting all the things 
together we have 

and, finally, 



Vs 



dlnV 
dk 

dlnV 



T'T H 



k 



a (T H - T)T ' 
dlnV a 2 T'T H 



(3.22) 



(3.23) 



din A; dk a (Th — T)T 

This is a general result which holds for any T(a) and a(t). One can further simplify ( 3.2SQ 
using the following redefinition: 



o(t) 



One can check that the expression for the spectral tilt then reduces to 



Vs 



H 2 



H 2 + H 



1 dT T H 



XY 



(3.24) 



(3.25) 



TdM (T H -T) 

where X (Y) denotes the term in the first (second) set of brackets. Note that while the 
first term depends on how the scale factor evolves with time, the second term primarily 
depends on how the temperature changes with the expansion of the universe, T(j\f). 
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3.5 Spectral Tilt and Amplitude during Bounce 



In this section we evaluate the magnitude of the spectral tilt in our scenario keeping in 
mind the various consistency requirements of the mechanism and the fact that we have to 
reproduce the correct amplitude of fluctuation in CMB. We recall that, as in inflationary 
cosmology, the cosmological fluctuations on a fixed comoving scale k freeze out at the time 
ti(k) when their wavelength crosses the Hubble radius. This is similar to what happens 
in the case of cosmological inflation. However, whereas in inflationary cosmology physical 
scales are increasing exponentially while the Hubble radius stays approximately constant, 
and thus scales are "expelled" from the Hubble volume, in a bouncing scenario such as the 
one presented in this paper, scales exit the Hubble radius in spite of the fact that during 
the relevant time interval the universe (and hence the physical scales) expands very little. 
Instead, the Hubble radius sharply decreases to a microscopic value starting from infinity 
(at the bounce point). 

Since the scale of the bounce is governed by A ~ M* which is expected to be of the 
order of the string scale, approximately in a string time the Hubble radius decreases to 
string length from infinity. All the scales which we observe today exit the Hubble radius 
in the process 9 . 

Let us first compute the temperature dependent factor Y in the equation ( |3.25| ) for the 
spectral tilt. As we will soon see phenomenological constraints will tell us that the bounce 
has to occur just after the onset of the Hagedorn phase when the background evolution is 
still governed by radiation (massless modes) as explained in the appendix ^. The task of 
computing Y is then easy as the temperature is just given by 

T = T b e~ M « T H e~ N . (3.26) 

Hence, 

Y = Th » Th , (3.27) 

T H -T T H -T b ' V ' 

since all the relevant modes are effectively expelled out of the Hubble radius very near to 
the bounce point. 

For (2.7) one can easily compute the first term, X, in the expression of the spectral 



tilt ( p-25 ), obtaining 



X = H2 = tanh 2 (At/^2) ^ A 2 t 2 

H 2 + H tanh 2 (AV\/2) + sech 2 (At/V2) ~ 2 + 0(A 4 i 4 ) ' K ' ' 

Putting everything together, we have a rather simple expression for the spectral tilt eval- 
uated at a length scale given by the comoving wavenumber k (note that the tilt in general 
is scale-dependent) 

^2^2 T}{ A 2 £ 2 T}{ 

7,3 ~ 1 = ^Th-T, ~ 2^ AT (3 ' 29) 



9 even if A is smaller than M* as required in our scenario, the picture remains true. 
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and we recall that the physical wavelength of the mode exits the Hubble radius at a time 
denoted by U(k). 

We will soon find out that for the physical scales that we observe in the CMB, U{k) « 
A -1 . In this approximation ti(k) is determined by: 

A , Xt X 2 ti{k) 

k = Ha = a = —= sinh —= Ft — - — (3.30) 
V2 V2 2 

so that 

The comoving wave number corresponding to the Hubble radius today is given by 

k Q = H a , (3.32) 

and thus we have 

_ 2 £ (3,3) 

Since during the bounce phase the universe expands very little, one can easily estimate 
ao ~ Tq/Tjj using the usual scaling during radiation era. Then it is easy to check that (see 
also appendix |B|) 

h s -l|-10- 60 f^V^ (3.34) 



V A J AT 

Thus we observe that the "bounce factor" X is typically very small and thus, as long as 
the temperature dependent factor Y is not too large we will reproduce a scale-invariant 
spectrum. 

For illustration let us consider a specific example which is consistent with all the 
requirements, viz. (i) we get a scale invariant spectrum, (ii) amplitude of fluctuation 
~ 10 as required by CMB, (iii) bounce occurs near the transition, i.e. pi, ~ M 2 ~ 
MpX 2 , and finally (iv) the fluctuations in energy are dominated by the closed strings in 
the Hagedorn phase. As explained in the appendix, the last constraint is met as long as 10 

AT 

— < 10" 30 (3.35) 
Th 

We will consider the limiting case, when AT /Th ~ 10~ 30 . Now the amplitude of fluctuation 
is approximately given by 



4® ~ 10- 10 ~ IJL^M S ~W- 10 M P (3.36) 

In other words an M s ~ 10 _10 M p , reproduces the required amplitude of perturbations 
observed in CMB. Requirement (ii) now forces A ~ 10 _10 M S giving rise to a very small red 
spectral tilt 

\rj s - 1| ~ 10~ 10 (3.37) 

10 We should point out that this is not an unnatural requirement. During the Hagedorn phase the 
temperature for most part remains very close to the Hagedorn temperature. The difference between the 
maximal temperature and the Hagedorn temperature is indeed exponentially suppressed. 
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This situation can of course change once we include additional corrections to the en- 
tropy relation, S(£ , v) near the transition, or if the scales involved are different. In fact, a 
really interesting case emerges if one considers the string scale to be around (Tev) which 
is also interesting from the point of view of particle physics phenomenology. In this case, 
A ~ 10 -15 M S ~ 10 -30 M p , and one obtains a spectral tilt which is just becoming observa- 
tionally significant! 

Before ending this section let us point out why for the phenomenological success, we 
require the bounce to occur close to the transition. If AT/Th -C 10~ 30 , for instance, then 
in order to get the amplitude of spectrum right, one has consider an unacceptably low 
string scale. Typically, deep inside the Hagedorn phase AT/Th is indeed very small, it is 
suppressed by an exponential of an exponential, and thus is not phenomenologically viable. 



4. Conclusions and Discussion 



In this paper we have extended the non-singular bouncing cosmology construction obtained 
in [Q] to more general matter sources, including a gas of closed strings. We have shown 
that this model leads to a simple implementation of the string gas structure formation 
mechanism proposed in [|[, provided that the energy density at the bounce is sufficiently 
large to lead a Hagedorn phase of the string gas, and the three large spatial dimensions 
are compact. 

The higher derivative gravity model discovered in is distinguished by being free of 
ghosts. No new gravitational degrees of freedom appear compared to those of Einstein 
gravity. The dilaton can be consistently taken to be frozen at all times. Thus, it is 
reasonable to assume that the metric fluctuations can be described using the equations of 
general relativistic perturbation theory. 

In this case, all of the conditions to obtain a scale-invariant spectrum of scalar metric 
fluctuations detailed in 0, |8| are satisfied. The matter fluctuations in the Hagedorn phase 
around the bounce time, when the Hubble radius diverges, are characterized by a heat 



capacity Cy which scales as r with the size of the region [11, y, Via the relativistic 
Poisson equation, these matter fluctuations induce a scale-invariant spectrum of metric 
fluctuations which propagate to late times, being squeezed, but without modification of 
the spectral shape. The spectrum has a red tilt, like in the string gas cosmology scenario 
H ||] . However, the magnitude of the tilt - calculated in this paper - is very small. 

Assuming that the initial state is chosen very early in the contracting phase, when the 
matter density is smaller than the current density, and assuming that the spatial size is at 
least one Hubble volume at that time, we obtain a cosmology which resolves the horizon 
and entropy problems of standard cosmology 11 . Thermal equilibrium over large distances 
is established in the contracting phase, thus allowing us to use equilibrium thermodynamics 
of strings in the Hagedorn phase, even if the latter is of a duration short compared to the 



11 By the entropy problem we mean the problem of obtaining a universe sufficiently large at the current 
temperature to include our current Hubble volume (see also |2?j] for another proposal to solve the entropy 
problem without inflation, making use of non-trivial dynamics of the extra spatial dimensions which string 
theory predicts). 
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physical length at the bounce time of perturbations relevant to current cosmology. Thus, 
our model provides an explicit toy model which demonstrates, in the context of an explicit 
action for the background space-time, that all the criteria required to obtain a scale- 
invariant spectrum in the scenario of || can be realized, and thus successfully addressed 
the concerns raised in ||. 
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A. Exact and Approximate Bounce Solutions 

In the previous work ^ , the presence of a cosmological constant was required for the late- 
time consistency of the ansatz (2.7) for the solution. To see this consider the behavior at 
large times. It is clear that the scale factor tends to a de Sitter solution: 

1 



a(t) 



while one can check that all the higher curvature terms vanish: 



PiHD 
^00 



£3 

n=0 



00 



seen 



w 



0. 



(A.l) 



(A.2) 



As a result, we are just left with Einstein's theory of gravity. In the absence of the 
cosmological constant, at large times when the higher derivative corrections become small, 
we will instead make a transition to the usual radiation dominated era of Einstein gravity. 
For the cosine hyperbolic ansatz the modified Einstein tensor looks like 



G 



oo 



Lq + L2 secri ( Ai/v2 ) + L4sech 



2) +L 6 sech 6 (\t/V2 



C + C 2 X 2 t 2 + C 4 A 4 £ 4 + 0(A°t' 



where the L's are given in terms of the function T by 



io(A 2 ) = ^ 
2 

9A^ 
2 

2 



L 2 (X 2 ) 
U{X 2 ) 

L 6 (X 2 ) 



1 -r(A 2 ) 



4^A 2 
M 2 



2^L A 2 + -r'(A 2 )A 2 



3A2 
2 



M 2 ' 

[r(A 2 ) 



1 - r'(A 2 )A 2 



(A.3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 

(A.8) 
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and the C's can be read off from L's: 



Co 

C-2 



Lq + L2 + L4 + Lq , 
-i(L 2 + 2L 4 + 3L 6 

C 4 = 4 0&2 + 6L 4 + 15L 6 
lb 



(A.9) 
(A.10) 

(A.11) 



and now we can match them with the coefficients coming from the stress energy tensor. 

Near the bounce point (i = and a = 1 by convention), we can expand the function 
p(t) in a Taylor series: 



Pit) = p\ 



t=o 



+ 



dp 

lit 



t + 



t=o 



2 dt 2 



t 2 , ifR 

+ 6 rf^ 



t=0 



,3 , J_^P 

+ 24 dt* 



t 4 + . 



(A.12) 



t=o 



The terms with odd powers of t vanish by our assumption of a symmetric bounce with 
bounce point t = 0. Replacing time derivatives by derivatives with respect to the scale 
factor, we obtain 



Pit) ~ P\ a =l + 



X 2 t 2 dp 



4 da 



4+4 



+ 



o=l 



"96" 



dp 
da 



+ 3 



d 2 



0=1 



da 2 



a=l 



(A.13) 



where we have used the specific ansatz (£? 

Thus, by matching the C 's in ( |A.4| ) with the coefficients of the stress energy tensor 
( A. 13; ), we essentially have three equations which determine the two dynamical quantities, 
Pb, A along with giving us a constraint. One can check that the exact solutions that were 
found in [Q] in the presence of a cosmological constant and radiation is a subset of the 
approximate solutions described above. 

In the case when one can apply the ideal fluid approximation fl2.9| ), the equations 
simplify considerably and one obtains 

3A 2 (12^A 2 + r'A 2 ) _ 2A 2 [6r - 7 - 4r'A 2 ] _ 48A 2 [42r - 43 - 30r'A 2 + 24c A 2 ] 

8 ~ 1 + uj ~ (1 +w)(9cj + 11) ~ 

(A.14) 

determining X 2 ,pb and giving rise to an additional constraint. One can rewrite the above 
equations in a more illuminating form: 

XP_ 

M 2 



Pb 

M 2 ' 



K 



l+UJ 



- 



A2 ) 



3A 2 ^ 



A 2 
IT? 



and 



Pb 



3M 2 X 2 G(^ 



where we have defined dimensionless functions: 



(— 

\M 2 



= 6T - 7 - 4FA 



'\2 



(A.15) 
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e (w) s -( 12 i£ +A2r ') < A - 16 » 

and 

/ A 2 \ 942r - 955 - 676r'A 2 + 576^ ? (m) 

fC ( ——k ) = - 2 ZL- - 3 ) *( A.17 

\M 2 J - 61PA 2 + 96r - 112 g 

A few comments are now in order. Since we want to restrict ourselves only to matter 
sources which are "consistent" (for instance ghost-free), they have to satisfy the dominant 
energy condition. This, in particular, implies 



p(l) > and ^ 
da 



< -1<u<1 (A.18) 

0=1 



imposing crucial restrictions on when one can have such a bounce and therefore prevent 
the singularity. Specifically, it says that the dimensionless functions that we introduced, 
viz. J 7 , Q, JC, has to be positive at the solution point. 

One may also wonder why we need to keep track of terms up to 0(t 4 ). Naively, one 
would expect that matching terms up to 0(t 2 ) should be sufficient to capture the physics 
near the bounce. We should remember that, for consistency, we need to solve the full 
Einstein equations, i.e. G mn = T mn up to 0(t 2 ). In terms of the "Hubble equation" (|2.5| ) 
this is equivalent to looking at terms up to 0(t 4 ). 

B. Modeling the transition from Hagedorn to Radiation 

We have seen in our analysis that for phenomenological success of the model, the bounce 
point should occur when string scale energy density is reached, i.e. right around the 
time when the massive modes are being excited and we are entering (or just entered the 
Hagedorn phase). In fact this is also important for the consistency of our computations 
of the perturbation spectrum, as we now illustrate. In order to be able to do any kind of 
"local" analysis with the metric fluctuations it is imperative that we are able to define a 
local matter density function, p(x), which is sourcing the metric. Now this is possible if 
the following limit is well defined: 

p(x) = lim^, (B.l) 

where M r and V r are the mass/energy and volume respectively inside a sphere of radius 
r around the point x. It is clear that a non-trivial and a well-defined limit only exists iff 
M r ~ r 3 which is indeed the case for ordinary particles. However in the Hagedorn phase 
E ~ r 2 and therefore the density function is ill defined and accordingly one cannot proceed 
to obtain its Fourier transform or power spectrum for that matter. Thus in order for us 
to perform the usual density perturbation analysis it is imperative that for the distance 
scales that are relevant for CMB, the energy is still dominated by the contribution from 
the massless modes which would give rise to an r 3 term 12 . 



As long as our detectors are sufficiently coarse grained, we will then have a well defined limit (B.l) 
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Now, in order to check whether this is true one would need to compare the energy in 
the massive string modes (Hagedorn matter) given by ( |3,71 ) 13 



E 



massive 



c hag r 2 M s 3 <j In 



with that of radiation (massless modes) 

-Erad 



1 - T/T H 



CradTV 



+ log log terms 



(B.2) 



(B.3) 



inside spheres of the relevant scales. Here c ra d,Chag ~ 0(1) constants. Now, in order for 
the limit (|B.1| ) to exist, minimally we require 



Ce 



E, 



rad 



E 



massive 



T 



M s r 



hi 



T 



1-T/T H 



> 1 



(B.4) 



where ro is the physical scale at the time of the bounce corresponding to the Hubble radius 
today. Therefore as long as the logarithm in the denominator ~ 0(1), since at the time of 



the bounce T ~ M s , to satisfy (B.4) we simply require 

Ce ~ M s r > 1 . 



(B.5) 



It is easy to estimate this quantity. We know that during the bounce the universe 
expands very little. Therefore, making the simplifying assumption that most of the expan- 
sion of our universe occurred during the radiation era since string scale energy density, we 
have 

— =>• M s r ~ — l - 



H, 



-i 



-f =► M s r ~ ~ 10 M » 1 . (B.6) 

IV1 s 1 

Indeed, ( (B.4j ) is easily satisfied for the relevant scales that we are observing in CMB today 
and therefore our analysis of perturbation spectrum is well justified. As an aside, we notice 
that for all scales r > ro the energy is going to be dominated by radiation and in particular 
the evolution of our universe would also be governed by radiation. 

What now becomes crucial is to verify that it is still the massive string modes of the 
Hagedorn phase that dominate the energy fluctuations in a sphere of radius ro, and not the 
massless modes. The amplitudes of fluctuation is proportional to the respective specific 
heats (more precisely the heat capacities) and therefore in order for the power spectrum 
( 3.13|) to hold we need 



Cc 



a 



rad 



a 



< i 



(B.7) 



Actually, from (3.7) one finds Shag, total 



massive 



-^massive 



+ k(M s r) 3 . However, the second term does not 
contribute to the fluctuations (heat capacity), neither does it pose any problem for the existence of the 



limit (B.l), since it scales as volume. In fact one sees that already the full expression of the Hagedorn 
energy suggests a regime (i? m assive fe(M s r) 3 ) when the fluctuations come from E mliSSive , but the energy 
is still dominated by the term linear in volume. Perhaps, this term could even be interpreted as the 
contribution from the massless radiative modes at the transition point. Unfortunately, we do not have a 
clear understanding of the thermodynamics during this transition phase and hence we had to "introduce" 
radiation by hand to capture the physics. Hopefully future research will clarify the situation further. 
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Differentiating the expressions for the energies with respect to temperature for the respec- 
tive phases we have 

Cc = r M s {l - £-) . (B.8) 
Thus, in order to satisfy ( |B.7| ) one requires 

AT T on , , 

— = 1 - — < 10" 30 . (B.9) 

C. Including the Dilaton and Other Light Moduli Fields 

So far we have not discussed the possible role of moduli fields, such as the dilaton or 
radion that are present in any string theory compactifications. It is well-known that for 
consistency of late time physics (fifth force constraints, and bounds on variation of physical 
constants) one typically requires 14 that these moduli fields are stabilized at a scale which 
is at least higher than the scale of big bang nucleosynthesis, i.e. Mev. There are two 
different possibilities that are compatible with our scenario, (i) If the scale of stabilization 
is higher than the string scale, then the bounce never probes this regime and the moduli 
fields always remain frozen playing no role either in the background evolution or in the 
generation of fluctuations, (ii) A second possibility, that the moduli are stabilized at a scale 
much lower than the string scale where all the dynamics is taking place is also consistent 
with our scenario. Let us look at this more carefully. 

In the Einstein frame, the action for these scalars just looks like a free theory 

1 r + 



St = --j d 4 v=swr, (ci) 

where eft has been appropriately normalized. An important subtlety is that after the confor- 
mal transformation (i.e. in the Einstein frame) the moduli fields such as the dilaton couple 
to the massive string modes. Thus, if there is stringy Hagedorn matter (massive closed 
string modes) present, then such matter would typically source the dilaton. However, as 
we explained before, we are primarily interested in a scenario where the energy density is 
not yet dominated by Hagedorn matter (only the fluctuations are) and in fact its energy 
density is hugely suppressed as compared to radiation. Therefore, such a source term can 
be completely ignored. In this context we note that since radiation is conformal, it does 
not yield any coupling to the dilaton after the conformal transformation. Let us therefore 
investigate the evolution of the free theory (C.l). 



Once we include the dilaton, the generalized Einstein equation reads 

^00 = i(/W + /v) (C.2) 

p 



where 



(C.3) 



14 It is possible construct some string theory motivated scenarios where one could avoid these observational 
constraints and still have a rolling light moduli field playing the role of a quintessence matter jj5|, which 
has become attractive due to the recent observations of dark energy. 
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This has to be supplemented by the dilaton equation 

<t> + 3H<j> = 0. (C.4) 
One can solve these equations. For instance, the dilaton equation yields 

P(j> = y = AM a ~ 6 • ( c - 5 ) 
Using the specific bounce solution we find 

7T 



m =4> b + y ^2 [tanh(t 7A72)sech(tyA72) + 2 tan" 1 (e'v 7 ^) - -j . (C.6) 

It is interesting to note that during the course of the bounce 4> changes only by a finite 
amount: 

A0 ^ 0(oo) _ ^.oo) = ttJ^. (C.7) 
Thus, we observe that if 

(C8) 

the string coupling 

9s ~ e*/ M * (C.9) 

changes very little. Moreover, since the kinetic energy of the dilaton is small, it does not 
effect the evolution of the metric. In particular, it is completely consistent to set the value 
of to be a constant. 
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